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I Applications of a simple thermal model to ultra-relativistic heavy-ion 

OO ■ collisions are presented. We compute abundances of various hadrons, in- 

(N ■ eluding particles with strange quarks, the p_L-spectra, and the HBT radii 

for the pion. Surprising agreement is found, showing that the thermal 
approach can be used successfully to understand and describe the RHIC 
data. 
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2h ' 1. Introduction 

■ 

With the wide stream of new high-quality data flowing from RHIC, as 
well as with the continued efforts at SPS (for recent results see, e.g., |2|, 
there is a growing need for a simple description of the basic underlying 
I physics. Only then our understanding of the phenomena occurring in ultra- 

relativistic heavy- ion collisions can be pushed forward, and space made for 
potential new phenomena, hitherto unexplained within the basic picture. In 
this lecture we argue that most of the "soft" features of the data from RHIC 
(particle ratios, momentum spectra, HBT correlation radii) can be explained 
very efficiently within an embarrassingly simple model, which merges the 
thermal model |, § |, |^, |, |, |lO|, |l5|, 0, 0, |l8|, |0 



X 



with expansion, and incorporates in a complete way the resonances 
[p3| , 24, 25, Our description uses hadronic degrees of freedom 

and starts at freeze-out, i.e. at the point of the space-time evolution of the 
system where the hadrons cease to interact. Pertinent theoretical questions. 
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such as what had been happening before freeze-out, what had led to the 
strong expansion of the system, why is the thermal picture successful, what 
is the nature of hadronization, not to mention the notorious "was there 
quark-gluon plasma?", will not and cannot be addressed in this lecture. 
Nevertheless, we believe that our studies prepare ground for such questions. 



2. The thermal model 

Historically, the ideas of the thermal description of a hadronic system 
go back to the works of Koppe [|^, Fermi Landau |^, and Hage- 
dorn [32|. More recently, many groups have used these ideas in numerous 



papers in an effort to explain the data from various relativistic heavy-ion 
experiments, from SIS, through AGS and SPS, to RHIC. Along the way, 
the original picture has been occasionally supplied with extra features, such 
as the fugacities controlling deviations from chemical equilibrium [B^] , finite 



volume and Van der Waals corrections M, 34 , or the use of the canonical 



instead of the grand-canonical ensemble [35, 3£, 37]. 

The works of Heinz and collaborators [^] put forward the concept of 
two freeze-outs. As the system expands and cools, it first passes through 
the chemical freeze-out point at temperature Tchem- Later, the particles can 
only rescatter elastically, until these processes are switched off at a lower 
temperature Tkin. In an appealing way the distinct freeze-outs explained 
the need for a higher temperature to reproduce well the particle ratios, and 
a much lower temperature to describe the slopes of the momentum spectra. 



In our work ||Tj, 23 , |2^] we have shown that with the complete treatment 
of resonances, the distinction between the two freeze-outs is not needed, 
at least for RHIC, and one can achieve very good explanation of all "soft" 
features of data assuming one universal freeze-out, 

^chem — ^kin — T. (1) 

We have also dropped, with the Ockham razor at hand, all other additions 
to the most naive thermal approach. The dropped features may be recon- 
sidered later on, provided there is a well-established phenomenological need, 
or theoretical argumentation. 

The main ingredients of our model are as follows: 

• There is one freeze-out, as discussed above, at which all the hadrons 
occupy the available phase space according to the statistical distri- 
bution factors. The scenario with a single freeze-out is natural if 
the hadronization occurs in such conditions that neither elastic nor 
inelastic processes are effective. An example here is the picture of 
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the supercooled plasma of Ref. |39|. Moreover, the STAR collabo- 
ration 41] has presented an important argument in favor of very 



weak rescattering after the chemical freeze-out at RHIC, based on 
the observation of the K*{892)^ peaks in the pion-kaon correlations. 
This essentially shows that the expansion time between the chemical 
and thermal freeze-out is shorter than the life-time of the i^*(892)'^, 
i.e. ~ 4 fm/c. Additionally, the fact that the measured yields of 
K*{892)^ pO| , PI are reproduced very well within the thermal model 
|ll3| , |l4| hints to the scenario with a short time between the two pos- 
sible freeze-outs, as proposed in Ref. |^]. Thus, approximation (||) is 
reasonable. 



A crucial feature of our analysis is the complete treatment of the 
hadronic states, with all resonances from the Particle Data Table |42| 
included in the analysis of both the ratios and the momentum spectra 
(cf. Fig. 1^). Although the high- lying states are suppressed by the 
thermal factors, their number increases according to the Hagedorn 




Fig. 1. A schematic picture of the hadronic soup, formed in an ultra-relativistic 
heavy- ion collision at freeze-out. The resonances decay subsequently into stable 
particles. The inclusion of many resonances is crucial for the success of the thermal 
approach. 
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hypothesis ||3^, 45, 46|, such that theh net effect is important. For 
instance, only a quarter of the observed pions at RHIC comes from 
the "primordial" pions present at freeze-out, and three quarters are 
produced via resonance decays. All decays, two and three body, are 
implemented in cascades, with the branching ratios taken from the ta- 
bles. For the p^-spectra the resonances are also very important, since 
their decays increase the slope, as if the temperature were effectively 
lower [14|. It has also been found that the inclusion of resonances 



speeds up the cooling of the system in hydrodynamic calculations |43| 



• Whereas expansion of the system does not alter the particle ratios at 
midrapidity (provided the system is boost-invariant, which is a good 
approximation at RHIC, see Sec. it becomes absolutely essential 
for the p^-spectra. We model the expansion (transverse flow) and the 
size of the system with two parameters: the proper time at freeze-out, 
T, and the transverse size, /Omax- 

• The model has altogether four adjustable parameters: two thermal 
and two geometric, which possess clear physical interpretation. The 
two thermal parameters, temperature, T, and the baryon chemical 
potential, fiB, are fixed by the analysis of the ratios of the particle 
abundances [0]. The two geometric parameters are fixed with help 
of the p^-spectra. The invariant time r controls the overall normal- 
ization of the spectra, while the ratio Pmax/T directly influences their 
slopes. 

All data used in the present study are for Au-|-Au collisions at ^/snn = 
130 GeV. 



3. Particle ratios 

The density of the ith. hadron species is calculated from the ideal-gas 
expression 

rii = Qi J d^pfiip), 
fi{p) = (exp mip) - fiBBi - fisS., - fiih) /T] ± , (2) 

where gi is the spin degeneracy, Bi, Si, and li denote the baryon number, 
strangeness, and the third component of isospin, and Ei(p) = ^p"^ + mf. 
The quantities //g, and /i/ are the chemical potentials enforcing the 
appropriate conservation laws. We recall that Eq. (^) is used to calculate 
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used for the fit, and 



Table 1. Optimal thermal parameters, ratios ^j^y^^ 

further predicted ratios. The preliminary experimental numbers for K*{892) |52 
have changed and better agreement with the model followed. 



the "primordial" densities of stable hadrons as well as of all resonances at 
the freeze-out, which later on decay. The temperature, T, and the baryonic 
chemical potential, fiB, have been fitted with the method to the originally 
available experimental ratios of particles, listed in the second group of rows 
in Table 1. The /is and fj,i are determined with the conditions that the 
initial strangeness of the system is zero, and the ratio of the baryon number 
to the electric charge is the same as in the colliding nuclei. It turns out that 
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the role of ^u/ at RHIC is negligible. 

For boost-invariant systems the ratios of hadron multiplicities at midra- 
pidity, dN/dy\y=o, are related to the ratios of densities, Ui, since 



dNi/dy 



dNj /dy 



y=0 



(3) 



The first equality follows trivially from the assumed boost invariance, while 
second one reflects the factorization of the volume of the system (see Sec. 



dNj/dy 
dNj/dy 



may be used to fit the thermal 



Hence the midrapidity ratios, 

parameters of the model. 

Table 1 presents the result of the fit. In our procedure the ratios mea- 
sured by different groups enter separately in the definition of x^- Thus, the 
number of the used data points is n = 16. The obtained optimal value of 
T = 165 lb 7 MeV is, most interestingly, consistent with the value of the 
critical temperature for the deconfinement phase transition obtained from 



the QCD lattice simulations: Tc 
and Tc = 173 lb 8 MeV for two massless flavors 
9 MeV lower than 174 MeV of Ref. 
obtained in Ref. |59]. 



154 lb 8 MeV for three massless flavors 
We note that our T is 
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|1|], and 25 MeV lower than 190 MeV 
Nevertheless, the results of the three calculations 
are consistent within errors. We have also computed other characteristics 
of the freeze-out: the energy density, e = 0.5 GeV/fm^, the pressure, P = 
0.08 GeV/fm^, and the baryon density, ps = 0.02 fm~^. We note that the 
results for the i^*(892) mesons, off by 50% when compared to the early 
preliminary data |52], came within the error bars of the data corrected later 
|0|, m. The lower part of Table | contains our predictions for particles 
containing strange quarks. The agreement with the data, released later, is 
very good. In particular, the triply-strange is properly reproduced. 

To summarize this part, we stress the high quality of the fit in Table 1 
for all kinds of particles, including those carrying strange quarks. 



4. Expansion 

Obviously, much richer information on the hadron production is con- 
tained in the transverse-momentum spectra. Various collaborations at RHIC 
measure, with impressive accuracy, the particle spectra of different hadrons, 
dNi/ (27rpj_dp±dy) , at midrapidity and for various centrality bins (the latter 
may be mapped to different impact parameters, |^0|). Unlike the case of 
the ratios of Sec. |^ modeling of the spectra involves not only setting the 
thermal parameters, but also a suitable inclusion of the expansion. Clearly, 
hydrodynamic flow modifies the spectra via the Doppler effect. Thus, an 
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important ingredient of our model is the choice of the freeze-out hypersur- 
face (i.e. a three-dimensional volume in the four-dimensional space-time) 
and the four-velocity field at freeze-out. Many choices are possible here, 
with some hinted by the hydrodynamic calculations. Our choice has been 



made in the spirit of Ref. pi] , |62| , p^ , p^ , 65, 66, 67, 68 1, and is defined by 

the condition ^ 

— — — = const. (4) 

Later on we denote the constant in Eq. (^) simply by r. In order to make 
the transverse size. 



P=\Jrl+rl, (5) 

finite, we impose the condition p < Pmax- In addition, we assume that the 
four-velocity of the hydrodynamic expansion at freeze-out is proportional 
to the coordinate (Hubble-like expansion), 

*^,^,^,^V (6) 



r T \ t t t _ 

Such a form of the flow at freeze-out, as well as the fact that t and 
coordinates are not limited and appear in the boost-invariant combination 
in Eq. (^), imply that our model is boost-invariant. We have checked 
numerically that this approximation works very well for calculations in the 
central-rapidity region. 

In practical calculations it is convenient to introduce the following pa- 
rameterization |]6^]: 

t = r cosh a|| cosh a_L, = r sinhay cosha^, 

= T sinh a± cos (j), ry = t sinh a±_ sin (j), (7) 

where ay is the rapidity of the fluid element, Vz = r^jt = tanhay, and 
a± describes the transverse size, p = rsinha_|_. The transverse velocity is 
Vp = taiaha±/ cosh ay. The element of the hypersurface is defined as 

dx^ dx^ 

dS, = W7^^^ dayda^#, (8) 

where = t, x^ = r^, = Vy, x"^ = and e^ai3y is the Levi-Civita tensor. 
A straightforward calculation yields 

dT,^{x) = u^{x) sinh(a_L)cosh(a_|_) da^da\\d(j), (9) 



such that the four-vectors dS^ and turn out to be parallel. This feature 
is special for our choice (§J^, and in general does not hold. 



8 



A question comes to mind as to what extent the assumptions (|^,^) are 
reahstic from the point of view of hydrodynamics. As a results of a typical 
hydrodynamic calculation, the freeze-out hypersurface contains, in the p-t 
plane, a time-like and a space-like parts |6^. 54. 65, |66|, 67, The 

' ] . Our param- 
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latter one is plagued with conceptual problems 
eterization neglects the space-like part altogether, thus avoiding difficulties. 
The time-like part of the hypersurface has, in many hydrodynamic calcu- 
lations, the feature that the outer regions in the transverse direction freeze 
out earlier than the inner regions. Our choice as well as commonly 
used versions of the blast-wave model, where the freeze-out occurs at a con- 
stant value of i, do not share this feature. On the contrary, our Eqs. (^J6|) 
correspond to the so-called scaling solution p2| , [7^ ] of hydrodynamic 
equations, which is obtained in the case where the sound velocity in the 
medium is low. Naturally, the validity of the assumptions and their rele- 
vance for the results should be examined in a greater detail. In Ref. |23] we 
have checked that two different models of the expansion lead to very close 
predictions for the momentum spectra at RHIC. Other parameterizations 
may be also verified with the help of the formulas given below. The fact that 
parameterization (0,^ works impressively well (cf. Sec. P), and at the same 
time the conventional hydrodynamic calculations have serious problems in 
explaining the RHIC data, hints, in our opinion, for a revision of the part 
of the assumptions entering hydrodynamic calculations and for extensions 
[75, 7f:, 77 1 of the picture used up to now. 



5. Decays of resonances 

The decays of resonances present a technical complication in the formal- 
ism. The resonances are formed on the freeze-out hypersurface with a given 
four-velocity. In the local rest frame of the fluid element the momenta of the 
resonances have thermal distribution, however, their decay products have, 
obviously, a different (non-thermal) distribution, which reflects the distribu- 
tion of the resonance and the kinematics. Below, we describe in detail our 
method, which is exact and semi-analytic (final expressions involve simple 
numerical integration rather than involved Monte-Carlo simulations). 

Consider a sequence of the resonance decays of Fig. ^. The initial 
resonance decouples on the freeze-out hypersurface at the space-time coor- 
dinate XN-, and decays after time r^v, with an average time proportional to 
the life-time I/Ttv- [| Let us track a single decay product. It is formed at 
the point a^Tv-i, decays again after time TAr_i, and so on. At the end of 
the cascade a particle with label 1 is formed, which is being detected. The 



^ In this section the indices i label the position in the cascade, and not the particle 
species, as in Sec. 0. 
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initial resonance 




final pion 



X 



freeze-out hypersurface 



Fig. 2. The cascade of resonance decays. 



Lorentz-invariant phase-space density of the measured particles is 
ni {xi,pi) = 

C^^P2„,„ „^ . r ^-ToTo /■j4„. ,(4)f^^ , P2r2 



(10) 



1712 



Xl 



d^PN 



B{pN,PN~l) / dTN^Ne 



PN 



X J dT.^ {xn) p% S^'^^ (^xn + - ^N^i^ In [pn ■ u (xn)] 



We have generahzed here the formula from Ref. |78] where a single reso- 
nance decay, without cascades, is taken into account. Note that the inte- 
gration over xn-1 ■■■ X2 is unconstrained, while the integration over xm is 
constrained to the hypersurface S. The delta functions impose the condition 
that the particle of velocity Pn/mn travels the distance from Xn to Xn-i in 
time Tn- The function B[q, k) is the probability distribution for a resonance 
with momentum k to produce a particle with momentum g, namely 

B{q.k) ^ 



Anp* 



E* 



11 



where b denotes the branching ratio for the particular decay channel, and 
p* (E*) is the momentum (energy) of the emitted particle in the resonance's 

^ In this notation b includes also the ratio of the spin degeneracies of the two particles. 
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rest frame. We note that B{k,q) satisfies the normalization condition 



J^B{q,k)=b. (12) 



Integration over all space-time positions in Eq. ([1^) gives the formula for 
the momentum distribution 

^vi-K^ = f d'^xini{xi,pi) = (13) 



d^pi 

^'^ B {p2,pi) ... [ ^^^J^B {pN,pN-i) I dT.^{xN) Pn fN\pN -uixN)] 



J Epj^ 

which should be used in the general case of any E and u. |^ 

We are now going to prove the second equality in Eq. (^). Starting 
from Eq. (^) we find the multiplicity of particles of type 1 coming from 
the discussed chain decay, 

/[ d^pN 
dT,^ (xn) / — — p% In [pn ■ u {xn)] , (14) 
Pn 

with an obvious notation for the branching ratios. Since the last integral in 
Eq. ( p^ ) yields an expression proportional to u^{xn), and the distribution 
function of the resonance is thermal, we can rewrite Eq. (|lj) in the 
equivalent form 



Ni = bN^N-l ...h^i j dT^^{xN)u^{xN)nN[T{xN),lJ'B{xN)\ 

= bN^N-1 ...h^inN {T, hb) / dT,fj,{xN)u''{xN). (15) 



Eq. (|T5|) indicates that the volume factor at freeze-out, / dS^ (xn) u^{xn), 
factorizes if the thermodynamic conditions (temperature and chemical po- 
tentials) are constant on the freeze-out hypersurface. This observation leads 
directly to the general conclusion that, as long as we integrate (measure) the 
spectra in the full phase-space (or, for boost-invariant systems, at a given 
rapidity y), the ratios of the particle yields are not affected by the flow 
and can be calculated with help of the simple expressions valid for static 
systems. This completes the proof of the second equality in Eq. (0). 



^ Note that the dependence on the widths Fk has disappeared, reflecting the fact that 
for the momentum spectra it is not relevant when or where the resonances decay. It 
is important, however, when and where the particles decouple from each other, which 
is determined by the choice of E(a;jv) and u{xn). 
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An important simplification follows if the element of the freeze-out hy- 
persurface is proportional to the four-velocity. This is precisely the case 
considered in our model where (compare Eq. (^) 

dT,^{xN) = (n:{xN)u^{xN)- (16) 



Then 



X / '^^-^B {pN,PN-l)PN ■ u{xn) JnIpN ■ u{xn)] 
^Pn 



dH (xn) Pi ■ u (xn) fi [pi ■ u [xn)] , (17) 

where we have introduced 

Pi-i ■ u (xn) [pi_i • u (xn)] (18) 

— {pi,Pi-l)pi • u(xAr) fi [pi ■ U (xn)] ■ 

The meaning of Eq. (|T^) is that as we step down along the cascade, the 
momentum distribution of the decay product, fi-i, is obtained from the 
momentum distribution of the decaying particle, with a simple integral 
transform following from the kinematics. In the fluid local-rest-frame, most 
convenient in the numerical calculation, we have u^{xn) = (1,0,0,0), and 
the transformation dlSh reduces to the form |14| 



/-i(^) = ^^/ dkkh{k), (19) 
2Egp*q Jk.(q) 

where the limits of the integration are k±{q) = mji \E*q it p*Eq\ /m\. Equa- 



tion (19) is a relativistic generalization of the expression derived in Ref. 
The technical advantage of Eq. (p^) is that the cascade can be performed 
in the rest frame of the original particle, with spherical symmetry and one- 
dimensional integrations over momenta, ([T9|), while in the general case of 
Eq. ( p!3[ ) only cylindrical symmetry holds and two-dimensional integrations 
over momenta must be used. 

In the case of three-body decays we follow the same steps as above, 
with a modification arising from the fact that now different values of p* 
are kinematically possible. This introduces an additional integration in Eq. 
(|l9|) . The distribution of the allowed values of p* may be obtained from the 
phase-space integral 

A I ^^^5{mR-E,,-E,,-E,,)5^^Hvi+V2+Vz)\M\\ (20) 

J ^pl ^PZ 
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where pi,P2 and P3 are the momenta of the emitted particles, Ep^^Ep^ 
and Ep^ are the corresponding energies (ah measured in the rest frame of 
the decaying particle), M is the matrix element describing the three-body 
decay, and ^ is a normalization factor. For simplicity we assume, similarly 
to that M. can be approximated by a constant, i.e. only the phase- 
space effect is included. Operationally, the final expression for three-body 
decays is a folding of two-body decays over p* with a weight following from 
elementary considerations based on Eq. ([20|). 

Finally, for the case satisfying condition (p!^), the spectra are obtained 



from the expression analogous to the Cooper- Frye ||73| , 31 1 formula, 

dN 



d'^p_\_dy 



j p^dY.^ h{p-u), (21) 



but with the distribution /i which has collected the products of resonance 
decays. With parameterization ^ we can rewrite Eq. (21) in the form 



T / da\\ / sinha^d (sinha_L) / d^p ■ u fi{p ■ u) , 

J-oo Jo Jo 

(22) 



d^p^dy j-oo " Jo Jo 



where, explicitly, 

p ■ u = m^cosha||Cosha_L — p± cos^ sinhaj,. (23) 

We end this section with a pedagogical discussion of the role played by 
various effects included. Figure |^ shows the m^-spectrum of positive pi- 
ous obtained with thermal parameters of Table 1. The dotted line shows 
the spectrum of primordial pious without expansion. The dot-dashed line 
adds the resonance decays; they contribute about 75% of the total, with the 
low momenta more populated. The dashed line is the result of the model 
with no transverse flow, i.e. including only the longitudinal Bjorken expan- 
sion. Finally, the solid line shows the full calculation, with resonance decays 
and the longitudinal plus transverse expansion produced by parameteriza- 
tion (§,^) . The characteristic convex shape is acquired as the result of the 
transverse flow. 



6. Transverse-momentum spectra 

Equipped with all elements of the model, we can now apply it to describe 
the p^-spectra. The thermal parameters are always those of Table 1. In 
principle, they could change with the centrality bin (impact parameter), 
but since the ratios of particles depend weakly on the centrality [|l|, |2|, |3| , so 
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Fig. 3. Contributions of various effects to the mj^-spectra of positive pions (normal- 
izations arbitrary, the relative norm of dotted and dash-dotted curves preserved). 



do the thermal parameters. We begin with presenting in Fig. 4 the fit to 
the ear hest- available minimum-bias data from the PHENIX collaboration 
[ |82| . We observe a very good agreement of our model with the data up to 
p_L ~ 2 or even, amusingly, 3 GeV. In that range the model curves cross 
virtually all data points within the error bars. At larger values of p±, where 
hard processes are expected to contribute, the model falls below the data 
for p and p. Since the values of the strange and isospin chemical potentials 
are close to zero, the model predictions for tt~^ and tt~ , as well as for 
and K~ are virtually the same. The value of the baryon chemical potential 
of 41 MeV splits the p and p spectra. Note the convex shape of the pion 
spectra. The vr"^ and p curves in Fig. 4 cross at p_L ~ 2 GeV, and the 
and p at p_L ~ 1 GeV, exactly as in the experiment. The values of the fitted 
geometric parameters are shown in second column of Table 2. 

The next plot. Fig. 5, shows an analogous fit made separately for 3 
different centrality bins. The obtained values of the geometric parameters 
are compared in Table 2. Again, the agreement is satisfactory]^ 

Finally, in Fig. |6| we show our results for all up-to-now available spectra 



^ For non-central collisions the shape of the hypersurface and the four-velocity at freeze- 
out is expected to be deformed in the x — y plane. In fact, in the hydrodynamic 
approaches this is the result of the elliptic flow, causing the azimuthal asymmetry 
of the spectra. The effect can be incorporated by properly extending the parame- 
terization (|^^. However, the effect of departing from the cylindrical symmetry by 
the amount needed to describe the elliptic-flow coefficient, V2, is negligible for the 
p^-spectra integrated over the azimuthal angle, considered in this lecture 64]. 
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1 2 3 1 2 3 



[GeV] 

Fig. 4. The p_L-spectra of pions (solid line), kaons (dashed line) and protons or 
antiprotons (dashed-dotted line), as evaluated from our model, compared to the 
PHENIX preUminary data obtained from Fig. 1 of Ref . [|| . Later official PHENIX 
data of Ref. |Q agree with the data used here. Feeding from the weak decays is 
included. 

at y^SAw = 130 GeV for the most-central collisions, including the particles 
involving strangeness. The upper part of Fig. ^ displays the spectra of pions, 
kaons, antiprotons, used earlier to determine the geometric parameters (last 
column in Table 2), and the predicted spectra of the (f) and i^*(892)'' mesons. 

The predicted spectrum of the (j) mesons agrees well with the reported 
measurement p^], with the model curve crossing five out of the nine data 
points. The (p meson deserves a particular attention in relativistic heavy- 
ion collisions, since it serves as a very good "thermometer" of the system. 
This is because its interaction with the hadronic environment is negligible. 
Moreover, it does not receive any contribution from resonance decays, hence 
its spectrum reflects directly the distribution at freeze-out and the flow. 
Thus, the agreement of the model and the data for the case of 4> supports 
the idea of one universal freeze-out. 

The upper part of Fig. |^ also shows the averaged spectrum of i^*(892) 
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resonances, with the data from Ref. [40|. Once again we observe a good 
agreement between the model curve and the experimental points. As al- 
ready mentioned in Sec. ||, the successful description of both the yield 
and the spectrum of K*{892)^ mesons supports the concept of the thermal 
description of hadron production at RHIC, and brings evidence for small 
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Fig. 5. Model vs. experiment for the PHENIX data |g3[ at three different central- 
ity bins for pions, kaons, protons and antiprotons. The thermal parameters are 
unchanged, while the geometric parameters following from the fit are given in the 
figure. 
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interval between chemical and thermal freeze-outs, in support of Eq. (|T]). If 
the i^*(892)'' mesons decayed between the chemical and thermal freeze-out, 
the emitted pions and kaons would rescatter and the i^*(892)'' states could 




0.5 1 1.5 2 2.5 3 3.5 



Pl [GeV] 

Fig. 6. (a) The p^-spectra at midrapidity of tt^, K~ , p, cj), and if*(892)°, and (b) 
of hyperons A, S, and fi. The asterisks represent the data from the STAR, and 
other symbols from the PHENIX collaboration. All spectra are for the most central 
collisions |8|, !§, H, H, H, |om, |8|l. The STAR data for the n' , R- , (j), 
S's and fl^s are preliminary. The updated experimental points for the S's are 
lower by about a factor of 2 from those shown, and much better agreement with 
the model follows. All theoretical curves and the data are absolutely normalized. 
The data and the model calculation include full feeding from the weak decays. 



17 



not be seen in the pion-kaon correlations. In addition, if only a fraction of 
the K*{892)^ yield was reconstructed, it would not agree with the outcome 
of the thermal analysis which provides the particle yields at the chemical 
freeze-out. Thus, the expansion time between chemical and thermal freeze- 
out must be smaller than the K*{892)^ life-time, r = 4 fm/c j40[ |. 

The bottom part of Fig. ^ shows the predictions of the model for the 
spectra of hyperons. Again, in view of the fact that no extra parameters have 
been introduced here and no refitting has been performed, the agreement is 
impressive. We note that the preliminary |87| data for the H's used in the 



figure were subsequently updated [^. The following reduction of the data 
by about a factor of 2 results in a much better agreement with the model. 
The data accumulated at lower energies at SPS showed that the slope of the 
O hyperon was much steeper than for other particles |9^]. On the contrary, 
in the case of RHIC the model predictions for the Q are as good as for 
the other hadrons. Since the Q contains three strange quarks, it is most 
sensitive for modifications of the simple thermal model used here, e.g. the 
use of canonical instead of the grand-canonical ensemble. The agreement of 
Fig. |6| does not support the need for inclusion of these effects. 

The various values of the geometric parameter sets are collected in Ta- 
ble 2. We also show their ratio, as well as the maximum and average 
transverse-flow parameter, /3, given in our model by the equations 

V ' ' rmax 

and 

We note that the ratio Pmax/T, and consequently, /JJ"^^ and {(3±), practically 
do not depend on centrality. 

To summarize this section, we conclude that the successful and economic 
description of the spectra, as seen from Figs. ^, and ^, provides a strong 
support for the thermal approach with universal freeze-out in the description 
of the ultra-relativistic heavy-ion collisions at RHIC. 



7. Excluded-volume effects 

In the presented model the fitted values for the geometric parameters, r 
and /Omax! are low, of the order of the size of the colliding nuclei. This leads 
to two problems: 1) the values of the HBT radii, as shown in Sec. ^ would be 
too small compared to the experiment, and 2) there would be little time left 
for the system to develop large transverse flow. Both problems can be solved 
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Table 2. The fitted values of the geometric parameters for various centrality bins, 
their ratio, and the maximum and average transverse flow parameters, as given by 
Eqs. (11,11). 



with the inclusion of the excluded-volume (van der Waals) corrections. Such 
effects were realized to be important already in the previous studies of the 
particle multiplicities in ultra-relativistic heavy-ion collisions [Q, |9l| , 
where they led to a significant dilution of system. In the case of the classical 
Boltzmann statistics, which is a very good approximation for our system 



[15 1, the excluded volume corrections bring in a factor [91| 

l + j:jVje-P^o/Tnj' 



(26) 



into the phase-space integrals, where P denotes the pressure, Vi = A^irrf 



is the excluded volume for the particle of species i, ^ and rii is the density 
of particles of species i. The pressure must be calculated self-consistently 
from the equation 

P = J2p^{T,^,, - Pv,/T) = 5:i^0(T,^,)e-^^»/^, (27) 

i i 

where P° is the partial pressure of the ideal gas of hadrons of species i. 
For the simplest case where the excluded volumes for all particles are equal, 



r, = r, Vi = V, the excluded-volume correction (26) produces a scale factor 
common to all particles, which we can denote by S^^. The formula ( p2|) 
becomes 

dNj n /■+°° rpra^yi/r r2-K 



(Ppj_dy 



T / da\\ / sinha_Lci (sinha^) / d^p-uS fi{p-u). 

J-oo Jo Jo 

(28) 



Tlie excluded volume per pair of particles is |7r(2ri)'^, fience tfie factor of 4 in the 
definition of Vi. 
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The presence of the factor in Eq. ( [2^ ) may be compensated by rescahng 
p and T by the factor S. That way, we retain ah the previously obtained 
results for the particle abundances and the momentum spectra. However, 
now the system is more dilute and larger in size. 

Next, we present an estimate of 5. With our values of the thermo- 
dynamic parameters we have J2i ^iC^^ f^i) — 80MeV/fm^, which leads to 
S = 1.3 with r = 0.6fm. Values of this order have been typically obtained 
in other works. Thus, the excluded- volume corrections can increase the size 
parameters at freeze-out by about 30%, and in consequence the problems 1) 
and 2) are alleviated: the geometric parameters become large enough to be 
reconciled with expansion, and the HBT correlation radii can be properly 
reproduced, see Sec. ^. 



8. HBT radii 

The transverse HBT radii -Rside and i?out (here we use the Bertsch- 



Pratt ||9^, Q parameterization) measured |^, Q at RHIC have 
values very close to those measured at smaller beam energies. Only the 
longitudinal radius, -Riongj exhibits a monotonic growth with ^snn p5| . 
The weak energy dependence of -Rside and i^out has come as a great puzzle, 
since the RHIC beam energy, y^J/viV = 130 GeV, is almost one order of 
magnitude larger than the SPS energy, y/SNN =17 GeV, and based on the 
hydrodynamic calculations one would expect much larger systems to be 
produced at RHIC. Also, a longer life-time of the firecylinder was expected 
at RHIC, which should be reflected in longer emission times of pions, which 
in turn would result in the ratio i?out/-Rside niuch larger than unity p8| . 
On the contrary, the experimental measurements indicate that i?out/-Rside is 
compatible with unity in the whole range of the studied transverse-momenta 
of the pion pair (0.2 < kx < 1-0 GeV). This fact is another surprise delivered 
by the analysis of the RHIC data for the pion-pion correlations. 

We have computed the pion HBT radii in our model. The calculation 



is based on the formalism of Ref. [99|, and is similar to the case of the 
particle spectra shown in Sec. Details will be presented elsewhere 100 |. 
The results of an approximate calculation neglecting the hadronic widths 
are shown in Fig. |^ where the HBT radii are plotted as a function of the 
transverse momentum of the pion pair, k^. We note that very reasonable 
agreement with the data is achieved. We have used S = 1.3 of Sec. |^. In 
particular, the ratio Rout/Rs\de (independent of the scale factor S) is close to 
unity. The k±_ dependence of -Rside is a bit too flat. The longitudinal radius, 
-Riong) is sensitive to the cut-off in the rapidity distribution and cannot be 
reliably computed in the present, boost-invariant, model. 
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9. Conclusions 



The presented results for the hadron production at RHIC support the 
idea that particles are produced thermally, and this is the basic lesson for 
today. The simple, economic model with the approximation of a universal 
freeze-out, simple expansion, and complete treatment of resonances, pre- 
dicts the particle ratios, the transverse-momentum spectra, and the HBT 
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Fig. 7. The HBT correlation radii for most-central collisions, i?side, ^out, and 
their ratio, as predicted by the model (solid line) and measured by the PHENIX 
collaboration. 
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correlation radii for the pion in agreement with the data. We note that 
the thermal approach works noticeably better at the RHIC energies than 
at lower energies, where, e.g., the particle ratios are not described to that 
accuracy |]l5|, or the spectrum of the $7 baryons is not reproduced. This 
indicates that the soft physics becomes simpler at RHIC, with our model 
being able to yield the quite impressive results of Table 1 and Figs. ^, ||, 
and |. 

In our phenomenology the pre-freeze-out stages are hidden and only the 
conditions at the moment where the hadrons decouple are relevant. This 
provides useful constraints for the more microscopic approaches. These cal- 
culations, describing early stages of the evolution, should ultimately provide 
the freeze-out conditions such as, or similar, to the ones used in our study. 

Certainly, the most challenging theoretical question which remains and 
should be addressed in future efforts is why the model works so nicely, and 
what it means for the underlying physics of particle production and the 
mechanism of hadronization. 

We are grateful Professor Andrzej Budzanowski for his encouragement 
and interest in this work, to Marek Gazdzicki for numerous helpful discus- 
sions, and to Boris Hippolyte for pointing out the early experimental spectra 
for the $7 baryons. 
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